Isotope effect in the pseudogap of high-temperature superconducting copper oxides 
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We study cuprates within Dynamical Cluster Approximation and find the pseudogap displays 
an isotope effect of the same sign as observed experimentally. Notwithstanding the non-phononic 
origin of the pseudogap the interplay between electronic repulsion and retarded phonon-mediated 
attraction gives rise to an isotope dependence of the antinodal spectra. Due to the strong momentum 
differentiation, such interplay is highly non-trivial and leads to the simultaneous presence of heavier 
quasiparticles along the nodal direction. We predict an isotope effect in electron-doped materials. 

PACS numbers: 74.72.-h, 71.10.Fd, 71.27.+a, 71.38.-k 



Copper oxide superconductors are characterized by 
the highest superconducting transition temperatures ever 
measured and, twenty-five years after their discovery the 
pairing mechanism has not yet been clarified. However, 
this is probably not the main reason why these materials 
are so intriguing. Indeed, due to the complex interplay 
between their electronic, magnetic and lattice degrees of 
freedom we still have a very incomplete understanding 
even of their normal (non-superconducting) phase. The 
most well-known anomaly characterizing the cuprates at 
small values of doping is the presence of a pseudogap. Us- 
ing the language of angular resolved photoemission spec- 
troscopy, the pseudogap is a strong suppression of the 
photoemission signal in selected parts of the Brillouin 
zone (BZ), occurring below a characteristic temperature 
called T* . This phenomenon is believed to be so relevant 
that any complete theory of cuprates should explain it. 

Differently from the superconducting transition tem- 
perature T c , T* displays a quite marked isotope effect 
characterized by a negative coefficient, i.e. it increases 
upon substituting Oxygen or Copper atoms with cor- 
respondingly heavier isotopes [l|-|4|]. This represents an 
unusual dependence, having conventional BCS theory in 
mind. Some possible scenarios for its explanation have 
been proposed: Assuming that the pairing in cuprates 
is directly mediated by phonons [H| so that the pseudo- 
gap is associated to the binding of two polarons with 
no long-range phase coherence, Ranninger calculated the 
isotope shift of T* and found good qualitative agree- 
ment with experiments. The main drawback of this 
approach is that a small electronic repulsion is enough 
to make bipolaron formation energetically highly unfa- 
vorable, and in cuprates local Coulomb repulsion can 
by no means be neglected. Another theoretical predic- 
tion available is based on the idea of a quantum critical 
point at optimal doping: By analyzing the effects of crit- 
ical fluctuations Andergassen et al. proposed a consis- 
tent picture of the experimental observation Q. Yet, no 
calculation has been done so far to see whether or not 
an isotope effect on T* is to be expected within proba- 
bly the simplest model for cuprates, namely a one-band 



tight-binding model for electrons or holes on the Copper 
sites with nearest-neighbor and diagonal hopping t and 
t' respectively, experiencing a local Hubbard repulsion U 
when occupying the same lattice site. Such model yields 
indeed one of the most successful descriptions of cuprates 
and many of the experimental observations such as d- 
wave superconductivity and the momentum differentia- 
tion characterizing the pseudogap have been qualitatively 
fairly well reproduced when solving it within Dynamical 
Cluster Approximation (DC A), Cluster Dynamical Mean 
Field Theory (CDMFT) and similar approaches @-[ll|- 

Yet, recent theoretical studies showed that in order to 
get a satisfactory quantitative agreement with photoe- 
mission and optics experiments the Hubbard (or the t-J) 
model is in fact not enough. Several experiments 



16( suggest a strong coupling to selected phonon modes 



and indeed things improve substantially if a coupling to 
phonons is included [17H23| . The size and the symme- 
try of such coupling have been previously calculated by 
combining ab-initio calculations and many-body effects 
24-281. Therefore we solve here the Hubbard model in 



the presence of electron-phonon (e-ph) coupling within 
DCA on 16-site clusters, as shown in Fig. [U We focus 
on its effects on the spectral function and look for the 
presence of an isotope effect in the pseudogap. 

Previous DCA studies of small clusters indicated that 
the pseudogap is very little affected and superconduc- 
tivity gets rapidly suppressed upon increasing the e-ph 
coupling 2j| 3^ . One may therefore expect our study to 
reveal hardly any isotope effect on the pseudogap This 
property in DCA is indeed a genuine electronic correla- 
tion effect and it is present already in the absence of 
phonons. On the other hand the different degrees of 
freedom are mutually entangled in cuprates so that even 
quantities which have no direct phonemic origin may dis- 
play an isotope dependence. Our study gives the novel 
evidence that this is indeed the case for the hole- and 
electron-doped (e-doped) cuprates: The Hubbard model 
has a pseudogap already in the absence of phonons and 
an e-ph coupling gives only a small change in its size and 
its momentum dependence. Nevertheless the pseudogap 



becomes isotope dependent consistently, as our results 
show, with experiments. The consequence of this is im- 
portant as the common way of looking at the isotope 
effect in cuprates is still based on the idea that either 
this is an extrinsic effect with little connection with the 
actual physical mechanisms or that it is a consequence of 
theories in which pairing is phonon-mediated, thus hardly 
realizable in cuprates. Our study shows that both of the 
above arguments are false: an isotope effect on T* is an 
intrinsic effect arising when a strong local Coulomb re- 
pulsion has to live together with a retarded attraction 
mediated by phonons and it is possible to describe the 
experimental data without necessarily invoking doubtful 
phononic pairing mechanisms [3l| . 
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FIG. 1: (Color online) (Left) 16-site clusters considered (the 
notation is the same as in Ref. [12]). (Right) Corresponding 
cluster momenta K c in the BZ. The K c marked by both the 
green full square and the red open circle are present in both 
geometries. The reddish and checkerboard gray areas denote 
the coarse graining regions for 16A and B respectively. 

The first question is therefore "What is the effect of 
phonons in the different sectors of the Brillouin zone?" . 
The form of the e-ph coupling term we consider is de 



rived from first-principle arguments following Ref. [25 . 
We show calculations for phonons of Holstein-type but 
we also considered breathing-type of modes in which the 
Oxygen atoms move in- and outwards the central Cop- 
per. For the DCA solver we have used is a Hirsch-Fye 
Quantum Monte-Carlo 



:-'».-, 



In Fig. [2] we show the spectral function for different 
cluster momenta K c for 5% hole- and e-doping. For hole- 
doping the pseudogap is in the K c = (it, 0) sector (panel 
a) while for e-doping the pseudogap is much smaller, 
though still somewhat visible at K c = (tt/2,0) (panel 
d) as well as at K c = (7r/4, 7r/2) 36|-Iq| . It is clear how 
in both cases the changes to the pseudogap features in- 
duced by the e-ph coupling are rather small. We observed 
that quantitatively such changes are fairly sensitive to the 
value of t' used, the cluster size as well as whether the 
e-ph coupling is increased beyond to the values shown 



FIG. 2: (Color online) 16-site cluster spectra at T=580K for 
different cluster momenta K c for hole- (a-b) and e-doping (c- 
d). l7=8|*|, |i|=400meV, \t'\=120meV. tv ph =200meV and 
lOOmeV for hole- and e-doped respectively. K c -resolved DCA 
spectra do not depend on any interpolation scheme. For the 
analytic calculation we have used Maximum Entropy [33Ll34|. 



here. Our definition of the e-ph coupling A is the ratio of 
the bipolaronic energy to the electronic bandwidth given 
by W=8|t|=3.2eV and we consider U—W. As we will see 
the point is that even if the A-dependence of the pseudo- 
gap is small, this is enough to give an appreciable isotope 
effect (i.e. a dependence of the pseudogap on the ion's 
mass, in other words on u> p h)- 

It is interesting to analyze the effect of the e-ph cou- 
pling on quasiparticles. For hole-doping these are in the 
K c =(7r/2, 7r/2) sector (Fig. [2)3) while for e-doping these 
are at K c = (-7r,0) (Fig. [2fc). At the temperature con- 
sidered here (580K) DCA gives in the hole-doped region 
much broader quasiparticles than for e-doping. The nat- 
ural consequence of this is that we can observe effects of 
the e-ph coupling in the e-doped side much better. In- 
deed in Fig. [2fc we see a substantial suppression of spec- 
tral weight (about a factor 2) while quasiparticles of panel 
b (hole-doping) are quite insensitive to A. The insensi- 
tiveness of hole-doped quasiparticles to A is found also if 
we look at spectral functions without "coarse-graining" , 
i.e. without the average over the BZ patches character- 
izing DCA [§] (sketched in Fig Q] as transparent regions). 
This procedure makes quasiparicle sharper, as indicated 
(for A=0 only) by the arrow in Fig. Yet, this is not 
enough to make e-ph effects discernible on the hole-doped 
side. Lower temperature studies with continuous-time 
QMC solvers (see Ref. |4l| for a review) are underway to 
understand whether this is more than just a temperature 
effect. 

The second question is then "Docs the e-ph interaction 
lead to an isotope effect in the pseudogap consistent with 
experiments?" . The answer can be directly read off from 
Fig. [3] and 0J In Fig. [3] the spectral functions at 5% 
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FIG. 3: (Color online) Pseudogap for different values of uj p h 
for 5% hole- and e-doping, upper and lower panel respectively. 



hole- and e-doping is shown for the K c -sectors where the 
pseudogap is observed. The outcome of DCA is clear: 
the pseudogap increases upon decreasing u> p h, i.e. upon 
increasing the ions' mass. This is the same trend found 
in experiments on hole-doped cuprates in which T* in- 
creases upon substituting 16 with ls O [H-Q. 

In order to extract more quantitative information from 
our data we look at a quantity which, unlike the spectral 
function, is independent on the analytical continuation 
method. This is the value of the K c -resolved Green's 
function at imaginary time r=/3/2. Indeed, this quantity 
does not depend on Maxcnt and it gives an estimate of 
the spectral weight at K c around uj = 0, Ak c (w). The 
results for different values of Lo p h are shown by the big- 
ger symbols in Fig. |4] For both the 16A and B cluster 
geometry ^K c =(ir,o)(0) decreases as uj p h decreases, as we 
expected from Fig. [3] Even though we considered a 
pretty coarse grid in 0J p h, from the slope of these curves 
we can extract an estimate of the isotope coefficient on 
^Kc=(7r,o)(0) for 5% hole-doping. The value we get is 
about 0.1. As it must be, the sign is opposite to that of 
the isotope coefficient on the pseudogap amplitude, be- 
cause a larger ^K c =(7r.o)(0) indicates a less pronounced 
pseudogap. Our study also allows us to make a pre- 
diction for what should happen on the e-doped side for 
which, to our knowledge, no isotope effect experiments 
have been performed hitherto. We expect them to re- 
veal a smaller isotope effect in the pseudogap (if visible) 
compared to the hole-doped side, but a substantial e-ph 
effect on quasiparticles. This is different from what pre- 
dicted by the resonant pairing scenario in which the 
isotope coefficient is essentially doping independent. 

More than one physical argument can be given to inter- 
pret our DCA results: Firstly, the trend of an increase 
in the pseudogap upon decreasing u> p h can be qualita- 
tively seen in terms of a simple paramctrization of the 



FIG. 4: (Color online) ^Kc=(7r.o) (w=0) for 5% hole doping for 
both clusters. The dashed black lines show the same quantity 
obtained for A=0 and changing U (oj p h~W0, 150 and 200mel/ 
corresponds to U =3.00, 2.95 and 2.90eV respectively). 



self-energy in terms of a pseudogap contribution plus a 
standard e-ph second-order diagram. An effect which 
may also play a role is a ph-induced broadening of the 
higher-energy features in the spectrum which gets smaller 
upon decreasing uj p h thus increasing the pseudogap [42| . 
Furthermore, our findings may be related to a recently 
proposed composite-fermion theory in which a reduction 
of the degree of coherence of quasiparticles due a smaller 
ui p h would lead to an increase of the "cofermion" binding 
energy and consequently of the pseudogap [43| ■ Here we 
however want to discuss the following very simple inter- 
pretation: Phonons mediate an attraction between elec- 
trons with a characteristic frequency scale set by u> p h- 
This introduces a reduction of the degree of correlation 
which is intrinsically Wp/j-dcpcndent, a fact that can nat- 
urally give an isotope effect on quantities which depend 
on correlation rather than on phononic mechanisms 44 1 . 
Our DCA results suggest that this mechanism can de- 
scribe the isotope effect in the pseudogap of cuprates. 
This can be seen in Fig. |4] where we plot ^4K c =(7r,o)(0) 
obtained by just slightly rescaling U keeping A=0 (small 
symbols connected by dashed lines). The reduction ob- 
tained by rescaling U is fairly similar to that obtained by 
really changing uj p h- Of course in our 16-site DCA cluster 
it is possible to describe a much richer behavior beyond 
this simplified picture: we indeed observe the ph-induced 
reduction of U together with a sizable mass enhancement 
of coherent quasiparticle. In Fig. [2] we not only see the 
tendency of the Hubbard band features to get closer to- 
gether upon increasing A but we also see the quasiparticle 
weight suppression, describing the typical loss of coher- 
ence induced by phonons in a metallic system. In other 
words we have a slightly less correlated system in the 
pscudogapped sectors of the BZ due to the ph-induced 
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reduction of U but heavier quasiparticles in the coherent 
sectors due to more standard e-ph effects. For a com- 
plete explanation of the isotope effect on the pseudogap 
it will be necessary to understand which correlations are 
directly affected by the e-ph coupling and, among these, 
which are the most important ones for the pseudogap 
formation. This remains an open question. We have 
shown that in the presence of a str ong differentiation in 
momentum space, like in cuprates [10(, the effect of the 
interplay between instantaneous repulsion and retarded 
attraction is highly non-trivial and leads to new effects. 
Since our calculation takes unbiasedly both ingredients 
into account and describes the experimental findings for 
the isotope effect in the pseudogap we conclude that such 
effects are not only relevant for cuprates but also for other 
strongly correlated systems in which the dominant mech- 
anisms are of electronic origin but some selected phoiion 
modes are still active. 
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